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SECTION 1.0 
INTRODUCTION 


During the past decade, satellite geodesy has pro- 
gressed to the stage where achievable accuracies are fre- 
quently comparable to or better than those obtainable by 
more conventional terrestial surveying methods, including 
both rod and chain triangulation and geodimeter traverses. 
This progress has been due in part to the development of 
more accurate instrumentation (particularly lasers), but 
also to refinements in the model for the geopotential field. 
Great progress has already been made in the estimation of 
spherical harmonic representations of the geopotential 
field, using "conventional" ground tracking data alone. 

When the mass of altimeter data from GEOS-3 is optimally 
utilized, substantially more accurate models are expected. 

Although satellite geodesy is perhaps most appro- 
priate for tying together points on a global scale, the 
results should be consistent with the normal surveys on 
the various datums throughout the world. Transformation 
to local (generally continental) datums from center-of- 
mass systems derived using satellite techniques has been 
a common method of checking the accuracy of results. His- 
torically, there have been two problems associated with 
this progress. First, the scale of the satellite results, 
and particularly the estimated station heights, has been 
largely imposed by the value of the geocentric gravitational 
constant which had to be supplied from non-satellite esti- 
mation processes - such as from launches of inter-planetary 
probes. At least unti], the last several years, the un- 
certainty in this constant has been on the order of 1 ppm 
and a significant source of error in estimation of the 
heights of satellite tracking sites relative to the earth’s 


center- 0 f-niass . The second problem, not unrelated to the 
first, is that comparisons of satellite derived coordinates 
with local survey coordinates have generally included a 
scale factor between the two systems, in addition to origin 
and orientation parameters. The result has been the util- 
ization of enough parameters in the transformation that 
satisfactory agreement between local datum and satellite 
derived coordinates was achieved. In the satellite estima- 
tion process, the only parameter whose scale should be 
considered uncertain is GM . But, when one attempts to find 
a value of GM„ for which the satellite results and the local 
surveys will be consistent, the result is so low C~398600.0km / 
sec ) that no one believes the result. A resoliition of this 
paradox is now apparently at hand, and includes the acceptance 
of problems with some ground surveys. 

Satellite derived positions are now approaching the 
accuracy of conventional surveys on all levels. For several 
years, Geoceiver derived coordinates with stated relative 
accuracies of 1 m in each coordinate have been claimed 
from the use of Navy navigation satellites. These positions 
have, however, also been subject to the scale uncertainty, 
with 5 m C.O.M. accuracy estimates. 

In this paper, we consider the utilization of laser 
and C-Band tracking of GEOS- 3 for the estimation of a set 
of global station positions, with the dominant set of track- 
ing from the continental United States. This set of positions 
is based predominantly on one day arcs of GEOS-3 during a 
10 day concentrated tracking period between February 23, 1976 
and March 5, 1976. These arc lengths are sufficiently long 
that station latitudes could be estimated for all stations 
and longitudes could be estimated for all except one reference 
station. All station heights were estimated. 


2 


SECTION 2.0 
SOLUTION TECHNIQUE 


The data set utilized consisted of the following 

arcs: 

7 1 day arcs of C-Band, laser, and S-Band data 

18 single pass arcs of data from lasers at Goddard, 
Bermuda, Grand Turk, and Patrick AFB, Florida. 

All passes dncluded tracking by at least 3 lasers 
On two passes, there was tracking by all 4 lasers 

All range data in the one day arcs was weighted with a’s 
of 1 m. All laser data was used, and the C-Band data was 
sampled at 1 data point per 6 seconds down to an elevation 
cutoff of 10®. Except for a few passes, a separate bias 
was estimated for each C-Band radar in each ai'c. This treat 
ment was used because it is conveniently implemented in the 
GEODYN data reduction program (T. Martin, et.al., 1976). 

Day to day bias stability of a meter or more has been dem- 
onstrated for a number of C-Band radars (Krabill and C. 
Martin, 1978). 

The objective in this weighting scheme was to obtain 
laser positions at the 4 calibration area sites with very 
accurate baselines, orientation, relative heights, and 
center-of-mass heights. The baseline, and most of the 
relative height information, was obtained from the heavily 
weighted single pass laser arcs. Since, however, the short 
arcs have practically no orientation or absolute height 
information, these are provided by the long arcs. The laser 
positions should thus be useful for the computation of 
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both very precise short arcs, and for the estimation of 
center-of-mass orbits of one or more revolutions in length. 

The gravity model used for the final station estimation 
solution was the GEM-9 model (Lerch, et.al., 1977). Several 
of the preliminary GEM-9 solutions were utilized in the data 
evaluation process prior to the availability of GEM-9. 

Overall station heights are determined by the value 

of GMq used in the solution. For lack of anything better, 

the value recovered with the GEM- 9/10 geopotential models 

was used. This value, 398600.64 km /sec , is the same as 

that estimated from lunar laser ranging data (WiHianis, 

3 2 

1974) and is within approximately 0.1 km /sec of other 

recent GM estimates. In the GEM-9/10 solutions, GM is 
e e 

determined primarily by LAGEOS tracking data. 


SECTION 3.0 

STATION POSITION RESULTS 


The estimated set of geodetic coordinates are listed 
in Table 1. All coordinates are independently adjusted 
except for 


1. The longitude of STALAS, which is constrained 
to the GSFC'73 value CMarsh, et.al., 1973). 

2. The relative coordinates of 

a. The Bermuda laser and C^Band radar 

b. The two Wallops C-Band radars 

c. Stations 4742 and 4452 on the island of 
Kauai (Hawaii) 

Separations between these pairs of stations were 
constrained to the local survey values. 

It will be noted that stations in close proximity 
(less than 4 miles) have been constrained to adjust together. 
This results in positions with a higher degree of internal 
consistency then would be the case with independent adjust- 
ments. It does, however, reduce the sol of checks which can 
be made on the validity of the solution. 
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Table 1. Estimated C-Band and Laser Coordinates 

Geodetic Coordinates* 


Station 
Name Number 

Geodetic Latitude 
(Deg, Min, Sec) 

Geodetic Longitude 
CDeg, Min, Sec) 

Height 

(Meters) 

STALAS 

7063 

39 

1 

13.3843 

283 

10 

19.7510 

14.96 

BDILAl 

7067 

32 

21 

13.8003 

295 

20 

37.8985 

-26.53 

NBEROS 

4760 

32 

20 

52.6323 

295 

20 

47.3765 

-18.88 

GRTLAS 

7068 

21 

27 

37.8189 

288 

52 

4.9867 

-22.10 

RAMLAS 

7069 

28 

13 

40.6886 

279 

23 

30.3059 

-27.05 

NWALI8 

4840 

37 

SO 

28.8957 

284 

30 

53.5420 

-29.98 

NWALIS 

4860 

37 

51 

37.0117 

284 

29 

26.4000 

-27.43 

BTR313 

4013 

26 

36 

56.2052 

281 

39 

7.7066 

-21.20 

ETRANT 

4061 

17 

8 

37.1790 

298 

12 

26.8721 

-5.07 

ETRMRT 

4082 

28 

25 

28.9359 

279 

20 

7.9908 

-23.20 

WSG219 

4150 

38 

58 

43.9386 

249 

S3 

20.1810 

1296.11 

WSM350 

4160 

32 

21 

24.3204 

253 

37 

14.9930 

1194.12 

WTRPPQ 

4260 

37 

29 

52.1543 

237 

30 

0.9544 

12.64 

WTRVAN 

4280 

34 

39 

57.0649 

239 

25 

6.8838 ' 

80.84 

WTRKPT 

4282 

21 

34 

19.7500 

201 

44 

0.1158 

298.58 

PMRPM4 

4446 

34 

7 

22.4154 

240 

SO 

42.9324 

-22.56 

PMRMR3 

4452 

22 

7 

57.9513 

200 

16 

25.5264 

503.85 

NELHAR 

4610 

39 

18 

30.5537 

244 

54 

48.0277 

2800.74 

WTRKAU 

4742 

22 

7 

24.7030 

200 

20 

3.7350 

1168.96 

KWAJM4 

4959 

8 

43 

17.6057 

167 

43 

36.0208 

4.43 

KALCOR 

4966 

9 

23 

55.2274 

167 

28 

55,9490 

59.84 

FRGMIO 

4960 

49 

8 

42.4363 

12 

52 

39.5260 

654. 55 

WSM398 

4198 

32 

25 

21.1072 

253 

36 

35.0614 

1183.94 


*a^ = 6378145 m, l/£ = 298.255 
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SECTION 4.0 

COMPARISONS ON THE NAD 


Several co'iiparisons were made between the estimated 
set of positions and local coordinates on the NAD 27. The 
first of those, shown in Table 2, compares the baselines 
between stations based on NAD coordinates, and based on the 
estimated set of center-o£-mass coordinates given in Table 
1. Only continental U.S. stations are listed. In addition, 
several continental U.S. stations have been omitted. The 
Merritt Island C-Band radar was not used because of its 
proximity to the Patrick laser, although the two stations 
were independently estimated. Similarly, only one of the 
White Sands main base stations was used. And the Ely, 

Nevada C-Band radar (Station 4610) was not used because of 
rather drastic revision of its NAD coordinates within the 
last several years, resulting in some question as to their 
validity. 

Based on the baseline comparisons in Table 2, it is 
evident that at least one of the station location sets is 
considerably in error. All but two of the estimated base- 
lines are larger than the NAD baselines, and a scale factor 
of 1.4 ppm is required to obtain the best fit transformation 
between the two datums. However, certain of the differences 
are known to be due to errors in NAD coordinates. The laser 
baselines will be discussed below. However, a geoidimeter 
traverse (Klosko and Krabill, 1974) comparison is available for 
the Merritt Island to Wallops baseline, with the results shown 
in Table 3. The precise traverse and C-Blnd estimated cords 
are in virtually perfect agreement. Such close agreement 
is no doubt accidental, but it does indicate that the East 
Coast distance discrepancies in Table 2 are largely due to 
NAD survey error. Also shown in Table 3 are comparisons of 


7 


Table 2. Comparison of Intersite Distances 
from Estimated Station Positions "w^ith Intersite Distances 
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Table 3. CoinparisoTi of Estimated Chord Distances 
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the Merritt Island - Patrick and Radar 398 - Radar 350 (at 
Wiite Sands) chords as computed from the NA.D survey and from 
the center-of -mass positions in Table 1. The agreements 
to within 55 cm at Patrick and 22 cm at White Sands are 
quite good, considering that these baslines were estimated 
without having any passes in the solution with simultaneous 
track on both ends of either baseline. 

The baseline distances in Table 2 between the East 
Coast stations and the IVhite Sands stations C4150 and 4160) 
are conceivably within the accuracy of tlie recovered solution 
The discrepancy in the station 4150 - station 4160 separation 
however, cannot be explained. The estimated baseline should 
be accurate to at least the 1-2 m level. If so, this points 
to a problem with either the surveys or some transcription 
thereof. 

The West Coast stations have posed a somewhat diffi- 
cult problem. The only really satisfactory tie among the 
estimated positions is the Vandenberg - Pt. Mugu 20 cm 
agreement. The others are worse than expected, both between 
themselves and with the East Coast stations. As shown in 
Table 4, the survey agreement is very good in relative lat- 
itude, and in the Pt. Mugu - Vandenberg relative longitude. 
The Pillar Pt. longitude disagreement is almost 4 m and has 
not been explained. The relative heights are also in dis- 
agreement at the 4 m level, considerably higher than the 
expected error in the station adjustment results. It may be 
noted that the "NAD” positions used for the comparisons 
are from the "Vandenberg adjustment", completed in 1971 
(Hieb, 1975), which included a geodimeter traverse along 
the West Coast. Agreement with the unadjusted positions 
would have been considerably worse. 
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In addition to possible survey problems and conceivable | 

problems in the recovery, the distortion between the NAD and f 

a center-of-mass system makes station comparisons in California | 

somewhat uncertain. Table 5 shows the computed shifts between I 

the estimated and NAD coordinates for the C-Band sites, and I 

comparable shifts for Geoceiver sites in the general area | 

(W. Strange, unpublished). These comparisons show a reasonably I 

smooth and consistent latitude shift, and an increasingly | 

negative longitude shift up the California coast. The longi- t 

tude shifts at Pillar Pt. relative to Pt. Mugu are approximately [ 

2 m larger for the estimated positions than for the Geoceiver [ 

positions. However, the Pillar Pt. longitude shift relative I 

to Vandenberg is about 4 m less for the estimated positions \ 

than the Geoceiver positions. The height shifts agree in | 

that Vandenberg should be shifted down the most. But the | 

recovered heights indicate that the shifts for Pt. Mugu and | 

Pillar Pt. should be approximately the same, whereas the J 

Geoceiver heights indicate that Pt. Mugu should be shifted [■ 

I 

down about 4 m more than Pillar Pt. This would make the I 

* i- 

height discrepancies in Table 4 even worse. 1 
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SECTION 5.0 
LASER STATIONS 


sSome comparisons are available for validating the 
relative positions of the laser sites. Dunn, et.al. (1978), 
have processed basically the same set of single pass laser 
arcs for GEOS- 3 that were used in our solution, with the 
objective of obtaining accurate baselines and relative 
heights. The baselines and relative heights from the Table 
1 solution are listed in Table 6, and the comparisons with 
the pure short arc solution is shown in Table 7. The agree- 
ment in baselines for the NASA lasers is well within the 
uncertainty of both solutions. Although a figure of 20 cm 
has been quoted, the la uncertainty is probably at the 10 cm 
level. Based on the sensitivity of the solution to changes 
in arc number, gravity model, etc., the relative height un- 
certainty is approximately 50 cm. 

The relative Patrick position is less well determined 
than the other stations, largely because only a few (6) arcs 
of data were available for use in the solution. For this 
station, we estimate the baseline uncertainty to be approx- 
imately 50 cm, and the la relative height uncertainty to be 
at about the same level. 



Table 6. Estimates of Baselines and Relative Heights 



















SECTION 6,0 

C-BAND STATION ACCURACY 


Because of the preponderance of simultaneous and near 
simultaneous track, the positions for Atlantic and continental 
U.S. stations are considerably more accurately determined 
than are those stations without simultaneous tracking by con- 
tinental stations. This is largely due to the minimization 
of orbit area in regions with large amounts of tracking and, 
by the same token, maximization of errors in regions of little 
tracking. Based on the sensitivity of the solution to the 
use of different gravity models, different methods of handling 
biases, different sets of data, and comparisons with other 
solutions, we estimate the accuracy of the adjusted stations 
to be 2 m for continental U.S. and Atlantic C-Band sites, 

5 m for the Hawaii sites, and 10 m for the German and Kwajalein 
radars. The laser relative positions are more accurate as 
discussed above. 

It will be noted that the solution has included a 
latitude adjustment for all stations, and that the comparison 
of the recovered latitudes with WGS 12 latitudes has shown 
approximately a 5 m difference. The GEM-10 (Lerch, et.al., 

1977) latitudes for the NASA lasers, however, show an average 
difference of less than 1.4 m from the recovered latitudes. 
Accordingly, the uncertainty figures above should be interpreted 
as including latitude as well as heights. Longitude uncer- 
tainty is relative to STALAS. 
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SECTION 7.0 
CONCLUSIONS 


Using one day arcs of concentrated C-Band tracking of 
GEOS- 3, station positions have been estimated which are be- 
lieved to be accurate at the approximately 2 m level and con 
siderably more accurate than most of the NAD surveys. The 
technique is capable of estimating all coordinates of all 
stations with the exception of a single reference longitude. 

The primary error source is geopotential model error 
and thus, as more accurate models are developed, accuracies 
(even from re-reduction of existing data) can be improved. 
C-Band results would certainly provide a good check on the 
results from the re-adjustment of the North American Datum. 



